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The local convertibility of quantum states, measured by the Renyi entropy, is concerned with 
whether or not a state can be transformed into another state, using only local operations and clas¬ 
sical communications. We found that in the one-dimensional Kitaev chain with quenched chemical 
potential p, the convertibility between the state for p and that for p + Sp, depends on the quantum 
phases of the system ( Sp is a perturbation). This is similar to the adiabatic case where the ground 
state is considered. Specifically, when the quenched system has edge modes and the subsystem size 
for the partition is much larger than the correlation length of the Majorana fermions which forms 
the edge modes, the quenched state is locally inconvertible. We give a physical interpretation for 
the result, based on analyzing the interactions between the two subsystems for various partitions. 
Our work should help to better understand the many-body phenomena in topological systems and 
also the entanglement properties in the Majorana fermionic quantum computation. 

PACS numbers: 71.10.Pm, 03.67.Mn, 05.70.Ln, 03.67.Lx 


I. INTRODUCTION 

Topological phases of quantum matter [1, 2], which 
cannot be described by local order parameters, have 
been an interesting subject in condensed matter physics 
because of the novel properties associated with these 
phases, such as the topologically protected ground state 
degeneracy [3], quantum anomalous Hall effect [4, 5], 
topological currents [6] and fractional quantum Hall ef¬ 
fect [7]. One prominent feature of topological phases is 
the possibility of creating nonlocal correlations among 
subsystems of the quantum matter. An ideal tool to 
characterize such correlations is the entanglement spec¬ 
trum (ES) [8]. It is a generalization of entanglement en¬ 
tropy and is defined as the eigenvalues of the entangle¬ 
ment Hamiltonian He, where He satisfies e~ HE = p (pis 
the reduced state of the subsystem) [9, 10]. Not only can 
ES be used to classify topological phases [11] but it can 
also detect the zero-energy edge mode [12, 13]. However, 
recent studies showed [14] that ES may not be universal 
for characterizing the quantum phases in the sense that 
it can exhibit singular changes within the same physical 
phase. Indeed, both a gapped system and a gapless one 
with distinct topological properties can have the same 
ground state [15], so that ES, when considered only for 
ground states, cannot distinguish between them. This 
also indicates that topological properties of a Hamilto¬ 
nian are not only related to the ground state but also 
the excited states and the energy spectrum. Therefore, 
it might be useful to consider ES that can reflect these 
factors for characterizing topological phases. 

On the other hand, in the quantum information com¬ 
munity, ES (or more precisely the eigenvalues of the re¬ 
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duced state of the subsystem) has been used to study the 
local convertibility of quantum many-body states. This 
topic is concerned with whether or not a quantum state 
can be transformed into another state, using only local 
operations and classical communications (LOCC). Let us 
concentrate on the most investigated bipartite pure state 
| i/jab) of a system divided into two subsystems A and 
B. As we know, under LOCC an entangled state can 
only be transformed into a state with the same or lower 
entanglement quantified e.g. by the von Neumann en¬ 
tanglement entropy [16] S v (pa ) = —tr(pAlogpA), where 
PA = trs(|i/)As)(V , AS |) is the reduced state of the subsys¬ 
tem A. However, the transformed state is not arbitrary. 
The question is: which states can be attained? This was 
solved for bipartite pure states through considering the 
Renyi entropy [17, 18] 

S a (pA) = rr— log(Trp^), (1) 

1 — a 

where a > 0. Consider two bipartite states | i/jab) and 
| iP'ab) w ith the corresponding reduced states pa and p' A 
respectively. If and only if S a (pA) — S a (p' A ) > 0 for 
all a > 0, then \iPab) can be transformed into | i/jab) 
through LOCC, possibly with the aid of a catalyst. Here 
the catalyst is an entangled state that participates in 
the transformation process but remains intact after the 
transformation. The Renyi entropy contains all the infor¬ 
mation of the eigenvalues of pa- For instance, Sq(pa) = 
log(d) (d is the effective rank of pa, be. the number 
of nonzero eigenvalues), lim Q _ ) .iS Q ,(pA) = Sv^a), and 
lim Q _ ) . oo S Q ,(/0A) = — logAi which is the single-copy entan¬ 
glement (Ai is the largest eigenvalue of pa) [19]. The 
Renyi entropy has been used to study the properties of 
the ground state | t/>(g)) of the XY model in the trans¬ 
verse field g [20]. It was shown that the differential 
local convertibility (DLC) of |V’(ff)) is related to quan¬ 
tum phase transitions. Here DLC refers to the transfor- 
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mation between | ip{g)) and ^{g + Sg)), where 5g is an 
adiabatic perturbation. In particular, for the transverse 
Ising model as a special case [21], DLC is affirmative in 
the paramagnetic phase (g > 1), while it can be negative 
in the ferromagnetic phase (0 < g < 1). It turns out 
that DLC is closely related to the size of the partitions 
and the correlation length of the system. When the cor¬ 
relation length increases to be comparable to the size of 
the subsystem, the Majorana fcmions (MFs), which form 
edge states, start recombining. As a result, DLC breaks 
down. 

However, for other models such as one-dimensional 
spin-i and spin-1 XXZ Hamiltonians, it was shown [22] 
that DLC is not necessarily related to quantum phase 
transitions but is a good detector of explicit symmetries 
of the model (e.g. the SU(2) symmetry of the Heisenberg 
model). In other words, DLC is more closely associated 
with the properties of the Hamiltonian than the ground 
state. This result is consistent with the conclusion that 
ES may not be universal for characterizing the quantum 
phases [14]. 

In this work, we study the local convertibility of the 
quantum states in the Kitaev chain [23] with a quantum 
quench. The quantum quench refers to the process that 
the Hamiltonian is changed abruptly so that the initial 
ground state is no longer an eigenstate of the new Hamil¬ 
tonian but evolves with time. Considerable efforts have 
been made to investigate the thermalization of quenched 
integrable and nonintegrable systems [24-37], as well as 
the dynamics of quenched topological edge states [38-45] . 
Quantum quench is essentially the natural dynamics of 
the system which was also studied for realizing quantum 
information processing in the spin systems [46] . The mo¬ 
tivation for studying the quantum quench in our work 
is that the state with time evolution involves many as¬ 
pects of the quenched Hamiltonian, including its excited 
states and the corresponding energy eigenvalues. Thus 
we conjecture that the local convertibility of such states 
may be in agreement with the quantum phases (espe¬ 
cially the topological properties) of the physical Hamil¬ 
tonian. Indeed, we found that our conjecture is true. 
The initial state for the quench is an uncorrelated state 
between the two subsystems. The property of DLC de¬ 
pends on whether or not the quenched system possesses 
edge modes, as well as the size of the subsystem. When 
the edge modes are present and the subsystem size for the 
partition is much larger than the correlation length of the 
MFs that form the edge modes, DLC is negative, other¬ 
wise DLC is affirmative. Our work points out the con¬ 
nection between the quantum quench and the topologi¬ 
cal properties of the system (the edge modes and MFs). 
Moreover, as the quench process simulates the quantum 
gate operations in the Majorana fermionic quantum com¬ 
putation (MFQC) [47], further investigation into DLC for 
the quench helps to better understand the entanglement 
properties of the quantum states in MFQC. 

The paper is organized as follows. The Kitaev model is 
introduced in Sec. II, where the numerical results regard¬ 


ing DLC is also presented. In Sec. Ill, we interpret the 
results of Sec. II, based on partitioning the system into 
two subsystems and analyzing the interactions between 
them. Finally, the conclusion is given in Sec. IV. 

II. THE MODEL 

We consider the one-dimensional Kitaev chain with the 
Hamiltonian [23] 

H = £ - w(cl n Cm +1 + h.c.) + A(c m c m +i + h.c .) 

m 

— h{ c ln C m ~~ 2 ); (2) 

where c' m , c m are the creation and annihilation operators 
of the electron in the mth site, w and A are the nearest- 
neighbor hopping and pairing amplitudes, and g is the 
chemical potential. The Hamiltonian (2) with periodic 
boundary conditions can be written in momentum space 
using Fourier transformation dfc = )Cm=i c m , 

f - 1 

= - £ (4, d- h )[R{k) • tr]{d k ,dl k ) T , (3) 
fc=- f 

where R{k) = (0, — A sin k, w cos k + g/2) is the pseudo- 
magnetic Held whose length times 2 is the one-particle en¬ 
ergy spectrum is 2 R{k) = yj (2 w cos k + g) 2 + 4A 2 sin 2 k, 
and a = ( tJ x ,<j y ,a z ) is the vector of Pauli matrices. 
Topological properties of the Hamiltonian was discussed 
in Ref. [44, 48] in terms of the winding of R(k) in the 
Ry-R z parameter space, where the region |/z| < 2 w is 
topologically distinct from \g\ > 2w in the sense that 
the winding of R(k) for the former surrounds the ori¬ 
gin of Ry-R z plane (A ^ 0), while it does not for the 
latter. It turns out that for |/z| < 2 w and A / 0, the 
Hamiltonian (2) with an infinite size and open boundary 
conditions supports a zero-energy edge mode composed 
of two unpaired Majorana fermions (a small energy split¬ 
ting for the edge mode is present for a finite-size chain). 
The Hamiltonian (2) is also equivalent to the long wave¬ 
length limit of the modified Dirac Hamiltonian [49] for 
describing the topological insulator with nontrivial Z 2 
index (=1). 

We consider a quantum quench problem as formulated 
below. The chain is initialized in the ground state |^q) 
for (w, A, g) = (w 0 , Ao, go). Then (w, A, g) are suddenly 
changed to {w',A',g') and the state will evolve accord¬ 
ing to the Schrodinger equation ihdt\^{t)) = 
i.e. l’I'(t)) = e~ ltH \ipo) (h = 1). The chain can be 
divided into two subsystems A and B , where A ( B ) 
comprises La {Lb) consecutive sites {La + Lb = N is 
the total number of the sites in the chain). The state 
l’P(t)) is in general an entangled state between A and 
B. If the initial state \ipo) is a product state between A 
and B, i.e. j^o) = | 4>a) <8> \<$>b }> then the entanglement 
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generated during the time evolution can be attributed 
solely to the action of the Hamiltonian (2). In this way, 
we envision that the entanglement properties of |H/(i)) 
will reflect some non-local properties of the Hamilto¬ 
nian. The initial product state, when restricted to be 
the ground state of the translation invariant chain (2), 
is either |0) 1 |0) 2 • ■ ■ 10)^ or |1) 1 |1} 2 • • • |1)^, where |0) m 
is the vacuum state of the mth site (c m |0) m = 0 and 
|l)m = c m|0)m)- The two states are the ground states 
for n —> —oo and p —> oo respectively. We also restrict 
ourselves to the situation that w = A and they are fixed 
during the quench. As the Hamiltonian can be scaled 
with H —>• H/w, we can essentially set w = A = 1. In 
summary, the scenario we consider is a quantum quench 
where p is suddenly changed from ±oo to some finite 
value / i' while w = A = 1 is fixed. The time evolution is 
written as 

=e- ltH ^' ) \vac), (4) 

where p’ indicates the chemical potential to be quenched, 
|vac) = |0) 1 |0) 2 • • • |0)jy is the vacuum (ground) state for 
p —>• — oo. In Appendix A, we shall show that the en¬ 
tanglement for the other case p —> oo is same to that for 
/i —>• —oo. 

We are interested in the asymptotic behavior of 
|\H(/z,t)) as t — > oo (p' is replaced by p for simplicity of 
notation). The question is: can the state |4/(/i + Sp, oo)) 
be attained from |tH(/x, oo)) through LOCC confined 
within the individual subsytems A and B (Sp is a per¬ 
turbation)? If the answer is positive, the quench process 
with the perturbation Sp can be simulated by LOCC and 
we say that | v h(/x, oo)) is locally convertible (or DLC is 
affirmative), otherwise the process involves non-local op¬ 
erations between A and B and we say that ^(/x, oo)) is 
locally inconvertible (or DLC is negative) [20] . The ques¬ 
tion can be solved through considering the Renyi entropy 
in Eq. (1). DLC is affirmative if and only if B^S^pa) > 0 
for all a > 0, or, 3^S a (pA) < 0 for all a > 0. Here, it 
is less important which state is convertible to which. We 
are more concerned with the overall convertibility. 

For the quadratic model (2) with quantum quench, the 
2 La eigenvalues of pa can be factorized as the tensor 
product of La vectors [50] 



( 5 ) 

where 1 > qi > q 2 > ■ ■ ■ > ql a >1/2 and they are the 
first La largest eigenvalues of the 2 La x 2 La correlation 
function matrix (CFM): C m<n = Tr pACmC' n with c m = 
( c m ,cl n ) T and 1 < m,n < La- The CFM is obtained 
from its Fourier transformation 

G(k, t) = [1 — R{k, t) ■ <r]/2 , (6) 

where R(k,t) is a time-dependent pseudo-magnetic 
field [44] associated with the quantum quench from R(k) 
to R'(k) defined in Eq. (3): R(k,t ) = cos(4 R't)R(k) + 



FIG. 1. (Color online). The sign of the derivative of the Renyi 
entropy with respect to the quenched chemical potential p for 
0 < p < 3 and 0 < a < 20. Here the unit of p is the hopping 
amplitude w which is set equal to the pairing amplitude A in 
the Hamiltonian (2). fa) The length of the subsystem A is 
L a = 10. (b) L a = 50. 

sin(4R't)fl(jfe) xR'(k) + [1 -cos(4 R’t)} [R(k)-R'{k)]R'{k), 

where R(k) = R(k)/R. and R!(k) = R'(k)/R'. When 
t —> oo, the oscillating terms in R(k, t) dephases and 

R(k,t) ->■ Reff(k) = -R' z {k)R'(k) with R'(k) = 
(0, - sin fc, cos k + p/2)/^TTJ^k + jRjZ, where the 
prime (') on the right-hand side of the equation is omit¬ 
ted for simplicity [51]. The parameters (w,A, p) for the 
initial state are technically chosen to be (0,0, —1). The 
Renyi entropy in Eq. (1) is simplified as 

1 La 

S a ( PA ) = -:-^log[g“ + (l-gj) Q ], (7) 

j = 1 

We numerically analyzed DLC through calculating the 
Renyi entropy for TV —»• oo and different La's. Fig. 1 
shows two examples with La = 10, 50. It can be seen 
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in Fig. 1(a) that for La = 10 and a fixed p £ (0,1.64], 
the sign of dfj,S a (pA) can vary with a, while for a fixed 
/i £ (1.64,3], dnS a (pA ) is always negative. This re¬ 
sult indicates that the quenched state is inconvertible 
for 0 < p < 1.64, but it is convertible for 1.64 < p < 3. 
Fig. 1(b) shows a similar result, but 1.64 is replaced by 
1.9. Namely, the range of p for negative DLC becomes 
wider. We also calculated other La's and found that 
indeed this range widens up to ( 0 , 2 ) as La increases. 
Further analysis of the data shows that for 0 < p < 2, 
DLC is affirmative when La is smaller than some critical 
value L c a as a function of p. When La > L A , the results 
are complex: for some ranges of p (e.g. 0 < p < 0 . 86 ), 
DLC varies with La up to a small number of incre¬ 
ments (Max(L^ — L c a ) = 6 ), but it is always negative 
as La increases further; for some other ranges of p (e.g. 
1.03 < p < 1.63), DLC is always negative as long as 
La > L a . For p > 2, DLC is affirmative for all La- 
The behavior of DLC is reminiscent of the edge mode of 
the subsystem A. The Hamiltonian of the subsystem A is 
Eq. (2) of La sites with open boundary conditions. When 
0 < p < 2, there are two Majorana fermions (MFs) 7 i, 72 
localized at the two boundary sites of the subsystem A 
respectively, with some overlap (a weak interaction H e ff) 
between them [23]. 

Heff = IT 7172 = e(c t c - i), rex e _iA/i °, (8) 

where <A = (71 — * 72 )/2, c = (71 + ij 2 )/2 are the creation 
and annihilation operators of the Dirac fermion (edge 
mode) formed by the two MFs, e = 2r is the energy of the 
Dirac fermion, and Lq = l/ln(^). When La Lq, H e ff 
is negligible and the two MFs are said to be unpaired, 
with the corresponding edge mode having zero energy 
(e —> 0). L 0 increases with p, and so does the overlap of 
the two MFs (for a fixed La)- When p approaches 2, Lq 
diverges (—>• 00 ) and the edge mode are absorbed into the 
bulk. When p > 2, no edge mode exists. Therefore, we 
can set a critical length [/ • Lq] such that an edge mode 
exists for La > \f • Lq] , where [x] denotes the minimum 
interger lager than x (La is an integer), and the factor 
/ designates the critical overlap of the MFs. The larger 
the factor / is, the smaller the critical overlap of the MFs 
becomes. In Fig. 2, L A and \f ■ Lq] are plotted against 
p with the optimized / = 2.582 for the best fitting of 
the two curves. It can be seen that the two curves fit 
well. In particular, the relative deviation 5L/\f ■ Lq] is 
smaller than 9% when p > 1.85 (see the inset in Fig. 2, 
SL = | L a — \f ■ Lq] I). The maximum deviation between 
L a and \f • Lq] for 0 < p < 1.99 is 7 sites. We shall 
discuss their connection in detail in the next section. 


III. PHYSICAL INTERPRETATION 

In this section, we shall interpret the result of the 
last section. To facilitate the discussion, we rewrite the 
Hamiltonian (2) in terms of the two subsystems A and 



FIG. 2. (Color online). (1) The main plot on a logarithmic 
scale: the solid pink line is L A versus p in the unit w(= A); 
the dotted blue line is \f-Lo] versus p, where [x] denotes the 
minimum interger lager than x, and / = 2.582 is optimized 
for the best fitting to the solid line and Lq = l/ln(2 /p). (2) 
The inset: the relative deviation 5L/\f ■ Lq] versus p, where 
SL = \L a — \f ■ Lq] |. 


B. 

H = Ha + H b + Hi , ( 9 ) 

where Ha (Hb) is the Hamiltonian ( 2 ) confined within 
the La (Lb) sites of the subsystem A (B). Hj = 
~w(c\ a c La+ 1 + cjyCi) + A(c La c L a+1 + cjvci) + h.c. de¬ 
scribes the interaction between the two subsystems (h.c. 
denotes the Hermitian conjugation of its previous terms). 
For w = A = 1, it is more elegant to write Eq. (9) in 
the basis of Majorana operators. Define c m = (a 2m -i + 
ia 2m )/ 2, cl = (a 2 m- 1 - ia 2m )/2 for 1 < m < TV, and 
let aj = /3j— 2 L A for 2 La + 1 5: J fs 2 (La A Lb) = 2TV 
to indicate that these MFs belongs to the subsystem B. 
We have 

La —1 La 

Ha — ^ ^ ^ ^ ^ ^2m-l^2m ? (10) 

m— 1 m= 1 

Lb — 1 Lb 

Hb — ^ ^ ^ /^2ra/^2m+l ^ ^ (^2m—l(^2m ? (H) 

m—1 m— 1 

Hi = ia 2 L A (3i + i/3 2 L B a i ■ ( 12 ) 

See Fig. 3(a) for illustration. Ha can be diagonalized as a 
sum of L a Dirac fermions [23]: Ha = Y)k =1 e k(^k^k — ^), 
where a k = (a 2k -i+ia 2k )/2, &j = w p m a m and W A 
block diagonalizes the coefficient matrix in Ha- Hb is 
similar: H B = Y)k~i e k (J>{h - jj), h = (fck-i+ifck)/^, 
Pj = Era W^mPm- 

We consider TV —> oo but La is finite (so Lb —>• oo). 
When 0 < p < 2, the subsystem B supports a zero-energy 
edge mode, say, ef B = 0, while for the subsystem A, 
only when La Lq does it support a zero-energy edge 
mode, say, e A —> 0 , as discussed in the context below 
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FIG. 3. (Color online), (a) The Kitaev chain with N sites is 
represented by a chain of 2N Majorana fermions. The vertical 
dashed lines mark the boundaries between the subsystems A 
and B. (b) When 0 < /z < 2, the subsystem A has an edge 
mode composed of two Majorana fermions (1 and 2) with a 
small overlap. The vertical dashed lines denote the modulus of 
the amplitude of the Majorana wavefunctions, and the dashed 
and solid curves are the overall shape of the wavefunctions. 
The subsystem B is similar. 

Eq. (8). In this situation, if we adopt the interaction 
picture, all the Dirac fermions oscillate except the zero- 
energy modes: dfc(f) = dk, bk(t) = e~ lUk bk- We 

expect that the main contribution to the quantum state 
in the infinite time limit will be the steady part of the 
interaction Hamiltonian, which has two types. (1) For 
the Dirac modes with non-zero energy, when ss e B , for 
some k, k ', Hj in Eq. (12) in the interaction picture will 
have a steady part oc a\bk> + h.c., which manifests the 
energy conservation. (2) For the zero-energy edge modes, 
the relevant steady part will be a linear combination of 
a,L A bk B , and their Hermitian conjugation. In 
terms of Majorana fermions, this part will be 

*X(«2L a ^ 2L b -1 + ^2L B ^2L A -l), (13) 

where x = l a 2 l a ' ^ 2 L b - 11 an ^ ^ le contributions 
relevant to W^ La _ 12La , W^ La} , W b Lb _ 12L? , W.? LbA is 
numerically found to be negligible and thus ignored. We 
shall argue that the second type of the steady part is 
the main source that renders the quenched state locally 
inconvertible. 



FIG. 4. (Color online). Entanglement (per site) of the 
quenched state for various length (La) of the subsystem A 
is plotted against ^ in the unit w(= A). Here the entangle¬ 
ment is measured by the von Neumann entropy of the reduced 
state of the subsystem A and ’’per site” means that the en¬ 
tanglement is divided by La- 

For Lb —> 00 and 0 < \l < 2, the edge mode of the 
subsystem B can be constructed analytically [23], 

02L b -i = sj 1 - Y ’ ( 14 ) 

3 = 1 

P2Lb = ^A ~ YX-^-'hLB-W • ( 15 ) 

3 = 1 

The energy of the corresponding Dirac mode = 
(/Ji + i$2L B )/^ is exactly zero. These results can be de¬ 
rived through examining the structure of the coefficient 
matrix W B . It can be seen from Eqs. (14,15) that the two 
MFs are localized at the boundaries of the subsystem B 
and their wavefunctions (the coefficient (—^ ) J_1 ) decay 
with j. The decay rate (= n/2) decreases with /z, indi¬ 
cating that the wavefunctions extend towards the inner 
part of the chain. The edge mode of the subsystem A, 
calculated numerically with finite La, has similar prop¬ 
erties. See Fig. 3(b) for illustration. With the extension 
of the wavefunctions, the number of the sites pertaining 
to the edge modes increases. We would expect that the 
effective dimension of the Hilbert space involved in the 
quench increases as well, and the entanglement becomes 
higher as the result of the larger dimension of the reduced 
state. However, this is not true, because the coefficient \ 
in Eq. (13), representing the interaction strength of the 
MFs, decreases (W b Lb _ 1 1 = y/1 — M 2 /4 from Eq. (14) 
and we found numerically that W^ La 2 l a is we ll approx¬ 
imated by yj 1 — /z 2 /4, so that x ~ 1 — /z 2 /4). 

From the above discussions, it can be seen that there 
are two competing factors that affect the entanglement 
properties of the quenched state when /i increases: ( 1 ) 
The wavefunctions of the MFs within one subsystem ex¬ 
tend inward; (2) The interaction between the MFs in 
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different subsystems decreases. Fig. 4 shows the entan¬ 
glement of the quenched state versus p for different La’s. 
At this juncture, we would like to remark that for a fixed 
p, the quenched state in the infinite time limit is also 
fixed, so that varying La amounts to changing the size 
of the partition for a fixed state. It can be seen in Fig. 4 
that the entanglement decreases very slowly for small p. 
When p is larger than some critical p c , the entanglement 
starts to decrease rapidly (e.g. in Fig. 4, p c ~ 1.6 for 
La = 10, and p c ~ 1.9 for La = 50). In terms of the 
above two factors, we expect that when p < p c , both fac¬ 
tors exist with the second one slightly more influential. 
When p c < p < 2, the wavefunctions of the two MFs 
in the subsystem A extend over the middle point of the 
subsystem and overlap considerably, so that the effective 
dimension of the Hilbert space involved in the quench 
saturates. Thus, the effect of the factor (1) on the entan¬ 
glement becomes negligible, while the factor ( 2 ) becomes 
dominant, causing a rapid decrease of the entanglement. 
When p> 2, the edge modes are absorbed into the bulk. 
It is conceivable that the entanglement decreases rapidly 
in this regime. This is because the bonding strength of 
the MFs within a single site of real space is stronger for 
larger p. See Fig. 3(a). As a result, the quenched state 
will be more localized to the single sites of real space, 
with lower entanglement. 

Comparing Fig. 4 with Fig. 1 and 2, we find that the 
value p c is consistent with the critical p for DLC to 
change. This indicates that for p < p c , although the 
entanglement, measured by the von Neumann entropy of 
the reduced state, decreases very slowly, the structure of 
the quantum state changes dramatically. The dramatic 
change is reflected in the fact that LOCC cannot simu¬ 
late the process in which p changes to /.i + Sp, that is, the 
quenched state is locally incovertible. Some insight into 
the inconvertibility can be gained by examining Eq. (5). 
Take La = 10 as an example. When the quenched pa¬ 
rameter p = 0.5, the sign of 9^ In qj roughly alternates 
between plus and minus when j increases (c^lngi > 0 ). 
This result influences the convertibility through the fol¬ 
lowing expression derived from Eq. (7) with log = In. 

La ^ 1 

d^S a (p A ) = ln( 9’ ( 16 ) 

3 =1 J 

where rj = (1 — qj)/qj and 0 < rj < 1. The sign of 
d^S a (pA) can vary with a (this is the case for p = 0.5 
here) on condition that the sign of 9 L In qj is not defi¬ 
nite when j varies. Namely, the uniformity of the sign 
of 9 M In qj is a sufficient but not necessary condition for 
the convertibility to hold. It can be verified that the for¬ 
mer is also a sufficient condition for majorization which 
is stricter than the convertibility [22], When p increases, 
the variation of the sign of 9 M In qj with j becomes less 
and less, while the convertibility remains to break down. 
When 1.64 < p < 2, the sign of 9 M In qj is almost uni¬ 
formly plus except a few of them (no more than 3 and 
j A 6 ). In this situation, the convertibility is restored 




FIG. 5. (Color online) The main plot: the ratio t c /x c versus 
La, where r c is the critical interaction between the two Majo- 
rana fermions in the subsystem A, as in Eq. (8), and Xc is the 
critical interaction between the Majorana fermion in the sub¬ 
system A and that in the subsystem B, as in Eq. (13). Here, 
“critical” means that when r/x > r c /xc, the quenched state 
will change from being locally inconvertible to being locally 
convertible. The inset: the ratio 7c/|AJj lid |“ versus La, where 
AJFj is the amplitude of the Majorana wavefunction in the 
middle point of the subsytem A when r/x = Tc/Xc- 


but majorization is not. When p > 2, the sign of 9 ^ lng^ 
is plus for all j, which guarantees the convertibility ac¬ 
cording to Eq. (16) and also the majorization. 

To justify the analysis, the ratio t c /xc is plotted 
against La in Fig. 5, where t c ,Xc are in Eqs. (8,13) 
with the subscript “c” added, and the critical p, L C A in 
the solid curve in Fig. 2 are used (the subscript “c” for 
r, x indicates “critical”). This ratio represents the com¬ 
petition between the above two factors, which can be 
seen from the definition of r and y. More explicitly, 
when t/x > t c /Xc, the quenched state will change from 
being locally inconvertible to being locally convertible, 
and we expect that T c /y c should be insusceptible to the 
change of La- Numerically, we find a slow decrease of 
Tc/Xc for large La, as in Fig. 5. For 100 < La < 500, 
T c /Xc ~ 0.08 ± 0.005. This is consistent with the data 
fitting in Fig. 2 as discussed below. As we know, r c is 
proportional to the overlap of the MFs. The overlap can 
be quantified as |A^ mid • A^ mid | = |A^ id j 2 , where A\ id , 
A c 2 mid are the amplitudes of the two Majorana wavefunc¬ 
tions respectively in the middle point of the subsytem A 
when t/x = Tc/Xc- Here A^ mid = A^ mid = A^ id . In the 
inset of Fig. 5, we find numerically that r c = l.llA^I 2 

for large La- Moreover, A(G d « yj 1 — (— i j-) i 5 '/ 2 and 

2 

Xc ~ 1 — x (see the context below Eq. (15)). Therefore, 
Tc/Xc ~ l-l(^) i A = l.le _z ^ /L °|^ Alc = l.le -2 582 = 
0.08, where L 0 is given in Eq. ( 8 ). 
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IV. CONCLUSION 

We have studied the local convertibility of the quantum 
state in the Kitaev chain with a quantum quench, where 
the chemical potential of the system is suddenly changed 
from ±oo to a finite value f.i. We found that the quenched 
state is locally inconvertible in the topological regime 
where the subsystems possess edge modes composed of 
Majorana fermions with weak interaction. When the in¬ 
teraction becomes sufficiently strong, or when no edge 
modes are present, the quenched state is locally convert¬ 
ible. However, the von Neumann entanglement entropy 
of the quenched state decreases for all the regimes, al¬ 
beit in the incovertible situation, the rate of decrease is 
very small. The distinguishing behavior of the convert¬ 
ibility of the quenched state, as compared with the en¬ 
tanglement entropy, indicates that the many-body quan¬ 
tum state may have rich structure that cannot be well 
characterized by the bipartite entanglement entropy. In 
particular, the rich structure, characterized by the lo¬ 
cal convertibility, turns out to be closely related to the 
topological properties of the system (the edge modes and 
Majorana fermions). Thus our work should help to bet¬ 
ter understand many-body phenomena in topological sys¬ 
tems, especially with a quantum quench. 

Future work can be pursued on more general param¬ 
eters, e.g. w ^ A in Eq. (2), in order to verify whether 
or not our conclusion regarding the edge modes and lo¬ 
cal convertibility is still applicable. Another interesting 
aspect is to examine a more general initial product state. 
For example, the initial state can be the tensor product 
of the ground states of the two subsystems. In the topo¬ 
logical regime where edge modes are present in the sub¬ 
systems, the quantum quench will involve interactions 
between the edge modes which amounts to performing 
quantum operations on the initial edge states. This is 
related to the Majorana fermionic quantum computa¬ 
tion where Majorana fermions are manipulated to realize 
quantum gates [47]. Detailed investigation will help to 
understand the entanglement properties of the quantum 
states in the Majorana fermionic quantum computation. 
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Appendix A: Identical Entanglement for quench 
from fi —» ±oo to ( w, A, ±/i) 

We shall prove that the quantum states with the fol¬ 
lowing four quenches have the same entanglement, i.e. 
the eigenvalues of the reduced states are identical. 

e -itH(w,A,±n) |f ul ), g—(ui,A,±/i) | vac ) ; (Al) 

where |vac) = | 0 > 1 | 0> 2 • • • 10 )^, |ful) = | 1 > 1 | 1> 2 ' ‘ ‘ | 1 >jv 
(| 0 ) m is the vacuum state of the mth site: c m | 0 ) m = 0 , 
and |l) m = <4|0) m ). 

Let us consider the equivalence between 
e -itH(wA,n)\f u \) and | vac ). The Fourier 

transformed correlation function matrix (CFM) for 
the former is G±(k,t) = [1 — Ri(k,t) ■ er]/ 2 , while for 
the latter it is G 2 (fc,t) = [1 — R. 2 (k,t) • er]/2. Here, 
Ri(k,t) = ~R 2 (k,t) = R(k,t ) whose detailed formula 
is presented in Ref. [44]. The eigenvalues of the reduced 
state of the subsystem A are the first La largest eigen¬ 
values of the CFM in real space C m , n = Trp^c m cf with 
c m = (c m ><4n) T - A s G\{k,t) and G 2 (fc,t) only differ 
on a minus sign, if the eigenvalue of G mj „ for Gi(fc,f) 
is (1 + A)/2, (0 < A < 1), there must be an eigenvalue 
(1 — A )/2 for G 2 (/c, t). Moreover, the eigenvalues of C m ,n 
can be written in pairs (1 ± A)/2 [50]. We conclude that 
the eigenvalues of C m ,n for G\(k 7 t) and G 2 (k 1 t) are 
same. This completes the proof. 

Next, we prove that e ~ ltH ( w -A,n) | V ac) has the same 
entanglement as e ~ ltH ( w A,-L) | V ac). As discussed above, 
the Fourier transformed CFM for the former is G 2 (fc, t) = 
[1 — i? 2 (fc,f) • cr ]/ 2 , while for the latter, it is G^(kA) = 
[1—i? 3 (/e, f)-<x]/2, where -R 3 (fc, t) corresponds to a quench 
process to (w, A, — fi). We notice that C m ,n for Ga(fc,f) 
is C mtn = f^ 7r G 3 (k,t)e lk( ' m ~ n ' > dk. If the integration 
variable k is shifted by 7 r, i.e. k —> k’ = k — n, we obtain 
an extra phase e l 7 d m_T d in the integrand. If this phase 
is negelected, the eigenvalues of C m ,n won’t change (this 
can be found by examining the entries of G mjrl ). There¬ 
fore, we essentially have G 3 (fc,f) —> G 3 (/c + 7 r ,t) = [1 — 
R 3 (k+n, t)-cr]/2. Furthermore, Rs(k+n,t) = —R 2 {k,t)- 
Thus, following the same argument in the previous para¬ 
graph, we conclude that the eigenvalues of the reduced 
states for e~ ltH ( w,A ’^\va,c) and e ~ ltH ( w A-v) | V ac) are 
same. The proof for the equivalence between the entan¬ 
glement of |ful) and that of e ~ itH ( w > A -^)|ful) 

is similar and omitted. 
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